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We consider two methods to find the effective parameters of the pion traversing a nuclear medium. 
One is the first order chiral perturbation theoretic evaluation of the pion pole contribution to the 
two-point function of the axial-vector current. The other is the exact, first order virial expansion of 
the pion self-energy. We find that, although the results of chiral perturbation theory are not valid at 

. . . normal nuclear density, those from the virial expansion may be reliable at such density. The latter 

CO ^ predicts both the mass-shift and the in-medium decay width of the pion to be small, of about a few 

O . MeV. 

O 

(N 

dc I. INTRODUCTION 

A considerable amount of work at finite temperature and chemical potential has been devoted to determining the 
effective parameters of strongly interacting particles in different media [1-6] . The results obtained are useful not only 
in analyzing the heavy ion collision experiments and properties of the early universe at different epochs, but also in 
extracting indications of an eventual phase transition. 
^ ' The case of pion appears to be the simplest to study. Being the Goldstone boson of the spontaneously broken chiral 
— ' symmetry of QCD, its interactions with itself and other particles are highly restricted by this symmetry, leading to 
' the effective theory of QCD, called chiral perturbation theory [xPT) [7,8]. At finite temperature, one has the further 
' advantage of having again only pions dominating the heat bath. Thus xPT provides a reliable method to calculate 
. the pion parameters at finite temperature [1-3]. 

' It is natural to apply xPT to calculate the pion parameters in nuclear medium [9,10]. Although the method parallels 
(<•-») , that followed for the case in a heat bath, the results calculated here to leading order may have restricted validity, due 
to the presence of baryonic resonances close to the -kN threshold. Similar difficulties also appear in determining the 
nucleon parameters at finite temperature [4]. 

In this work we compare the result with that of the (first order) virial expansion of the pion self-energy 

O ' [4,11,12]. The latter gives the shifted pole position in terms of an integral over the product of the density distribution 
^ . function times the TriV scattering amplitude obtainable from experiment. It is thus free from the difficulty encountered 
^ ' in calculating the amplitude and is valid as long as the nucleon 'gas' is dilute enough. 

The virial formula has also been employed earlier to the same problem, but only in an approximated version [13-15]. 
If, however, the amplitude varies appreciably in the range of integration, in particular, if it changes sign - as is the 
, case here - this version is not justified. 
■ We first rederive the results of xPT by evaluating the axial-vector current correlation function to one loop, using 
the in-medium Feynman rules for the original chiral Lagrangian in presence of external fields. Besides completeness, 
its purpose is to show that this conventional framework is quite simple, without requiring functional integration over 
the nucleon field to produce a 'new' effective Lagrangian [10]. We then derive the exact, first order virial expansion 
for the pion self-energy and evaluate it with experimental data. 

Sec. II reviews briefly xP^i constructing the Lagrangian for the ttTV system [16,17]. In sec. Ill we work out the 
shift in the mass and the decay constant of the pion using this Lagrangian. Next we derive the virial formula and 
evaluate the pole shift in sec. IV. In sec. V we discuss the limitations of these methods. 

II. CHIRAL PERTURBATION THEORY 

The Lagrangian of QCD with two massless quark flavors is 

C^QlD=^nrd,q + ---, ^-{i] > (2-1) 
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where the dots denote terms involving other fields. If we split the quark field into its right and left handed parts, 
<1R,L = 5(1 ±75)0', it is clear that C^qcd is invariant under the symmetry group G = SU{2)ii x SU{2) l of independent, 
global SU{2) transformations on and ql, 

qR^gRQR , qL gLQL , 9r,l & SU{'2)r^l ■ (2.2) 

Phenomenology suggests strongly that the symmetry of the Lagrangian is broken spontaneously by the vacuum state 
to the diagonal subgroup H ~ SU{2)v, giving rise to the pionic degrees of freedom. 

In xP'r 0116 derives the transformation rules for the observed Goldstone and non-Goldstone fields from the above 
symmetry of the underlying QCD theory. It turns out that the Goldstone fields n^{x), i = 1,2,3 are collected in the 
form of a unitary matrix 

3 

u{x) = e'^(^)/2^' , Tr{x) = ^7r'(a;)r' , (2.3) 

where the constant can be identified with the pion decay constant, = 92.4 MeV and r* are the Pauli matrices. 
Then the matrix u transforms under G according to 

u QRuh'^ = hugl , (2.4) 

where the group element /i(7r) G SU{2)v- Notice that h is .x-dependent due to its dependence on 7r*(a;). However the 
square of this matrix = U has the global transformation rule 

U^gnUgl (2.5) 

On the other hand, the non-Goldstone, nucleon doublet field tp{x) transforms as 

^^h^ , V'=(^) ■ (2.6) 

There are two Noether currents following from the symmetry of jCqqjj, namely the vector and axial vector currents, 

V;{x) = q{x)j^-q{x), A'^{x) = q{x)^^l5^q{x) . (2.7) 

The evaluation of the correlation functions of the currents is most conveniently carried out in the external field 

method [8] . Although we are interested here in such a function of the axial- vector current only, we couple both the 
currents to the external fields v'^^^{x) and a^(a;) to reveal the full symmetry of the underlying theory. Thus the original 
Lagrangian extends to 

C%lD+v],{x)Vt{x) + al{x)A>;{x) 
= iqRl^id^i - i{Vf, + a^)}gR + iqLl^{d^ - i{v^, - a^,)}qL H (2.8) 

where v^{x) and a^{x) are the matrix valued external vector and axial vector fields, 

3 i 3 j 

^^^(a;) = I^<(a;)y , a^(x) = ^a;(a:)y. (2.9) 

i=l i=l 

The extended Lagrangian is now invariant locally, i.e. under x-dependent symmetry transformations on qn and qL, if 
the external vector and axial-vector fields are also subjected to the appropriate gauge transformations. 

The presence of the external fields in the underlying theory and the associated gauge invariance can be readily incor- 
porated in the effective theory. All we have to do is to replace the ordinary derivatives by the covariant derivatives [18]. 
Thus we have for U , 

D^U = d^U - i{v^ + a^)U + iU{v^ - a^), (2.10) 

and for tp, 
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'i^ij. = -^(u^[di^-iivij. + a^)]u + u[dij,-i{vij,-a^)]u^^ . (2.11) 

The building elements for the effective Lagrangian at this stage are thus U, Df^U, tp and V^xj}. 

A simplification in the construction of the effective Lagrangian emerges by noting that the variables {U , DfJJ) 
transform under the full group G, while 'Dnij)) transform only under the unbroken subgroup H. We may take 
advantage of the mixed transformation property of u to redefine the former type of variables so as to transform under 
H only. Thus one introduces the variable [16], 

u^ = iv)D^Uu^ =ul , (2.12) 

replacing U and Dfj,U [18]. Any term in the Lagrangian that is constructed out of these variables so as to be invariant 

under H will also be automatically invariant under G. 

We now write the effective Lagrangian of X-f ^ for ttN system as 

^eff = + ^N, 

where is the well-known pion Lagrangian [8] , which to leading order is given by 

C^ = ^{{D^UD>^U^)+ml{U + U^)} , (2.13) 
(• • •) denoting trace over the 2x2 isospin matrices. The pieces in £jv to first and second order, 

= jC^^^ + (2-14) 

are 

= ^^p- ruN)^ + Y'PM , (2-15) 

and 

£j^^ = c^ml{U + U^)^iP - ^{u^u^){^V^'V''xlj + h.c.) + ^{u^u^^iP - "-^i^^YKMi^ ■ (2.16) 

We shall use vertices in to second order and those in to first order in our perturbative calculations. Here 
QA turns out to be the axial-vector coupling constant appearing in the neutron beta decay, (/a=1-27. The coupling 
constants Ci , C2 and C3 are determined from the experimental data for -kN scattering in the low energy region and its 
extrapolation inside the Mandelstam triangle, where it is compared with the xPT evaluation, getting [20,17], 

ci = -0.81 ± 0.12 GeV-\ C2 = 3.2 ± 0.25 GeV-\ C3 = -4.66 ± 0.36 GeV"^ . 

Expanding out in the pion field and setting v^=0, we bring out explicitly the vertices, contributing to the pion pole 
diagrams to one loop [19], 

= - F^d^T, . a'* , 

.2., _ , , 



(2) 2m^ci - C2/I. „ 2 y 

+C3 (^;^5^7r • S'^TT - -|-a^7r • a^^ • (2-17) 



where the superscript '0' indicates free Lagrangian densities. 
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III. MASS AND COUPLING SHIFTS 



The two point correlation function of the axial vector current in a medium is given by the ensemble average, 

i J d^x e'«-^ Tr [e-^^^-^'^^Alix) (0)]/Tr [e-/3(^-/^^)] . (3.1) 

Here H is the QCD Hamiltonian, /3 is the inverse temperature and A'^ is the number operator for nucleons with 
chemical potential fx. Our aim is to find the corrections to the pion pole in this correlation function due to interaction 
of pion with nucleons in medium at zero temperature. 

We shall work in the real time formulation of the field theory in medium [21]. Here the perturbation expansion 
proceeds as in the conventional (vacuum) field theory, except that the propagators assume the form of 2x2 matrices. 
For the pole term to one loop, it however suffices to work as in the vacuum field theory, only replacing the vacuum 
propagators by the 11-component of the corresponding propagators in medium. 

So we first calculate the vacuum correlation function 

i J dSe^«-- {0\TA\{x)AUom , (3.2) 
in xPT. We recall that the generating functional of QCD, 

(0|Te'.'' '*^^"M(^)^i'(^)|0) , (3.3) 

is represented in X-P^ by 

(OlTe'-'' '''''^'"'(''■'^■"^lO) , (3.4) 

where Cint is obtained from Eqs.(2.17). Since the two point function (3.2) is the coefficient of the term quadratic in 
a^(x) in the expansion of the generating functional (3.3), we find these quadratic terms from the functional (3.4) of 
the effective theory. 

(a) /WW WW 



(b) ^V\JW-(^-^AfVV\ (c)/WW--^— WW 



(d) WAA<^---WW (e) WW----^]5WW 



WW (9)/WW 

FIG. 1. Feynman diagrams for the two point function to one loop. Only loops with nucleons are considered. 

Fig. 1 shows the free pion pole diagram (a) giving the amplitude, 

S"qxqatF^A{q) , A{q) = l/{q^ - ml + le) , (3.5) 

together with all one loop corrections to it, relevant in the nuclear medium. Let us calculate the self-energy diagram 
(b) to illustrate the method. Its contribution may be written in the form 

6'''qxq.tF^A{q){-iU{q)}A{q) , (3.6) 

which modifies the free pion pole term (3.5) to 

q^ -ml- Il{q) 
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To calculate the self-energy function n(g) of the pion in the nuclear medium, we first write it in vacuum, 

n<°H9) = ^J + m^)i^^{J, - i + m^)] A^(p)A^(p - q) , (3.8) 

where A^(p) is the vacuum nucleon propagator after extracting the factor {^ + 171^), {p) = i/{p^ — nij^ + ie) . 
Following our discussion above, we now replace the vacuum nucleon propagator in Eq. (3.8) by its 11-component in 
nuclear medium {Ep = y^p ^ _|_ m|f), 

^n(P) = p2_^2^ + i, - 27T{n+iEp)e{po) + n-iEp)e{-po)} Sip^ - ml) (3.9) 
where n^{Ep) are the distribution functions respectively for the nucleon and the antinucleon, 

-^(^^^= e^iE^L + l - ^3.10) 
As the temperature goes to zero, we get for /i > 0, 

n+{Ep)-^9itx-Ep) , n-iEp)^0. (3.11) 
The density dependent part of the self-energy in the medium is then obtained as 

n id) - -Sg^jjq j -4(p . qr + 1' ■ ^ ^ ^ 

Consider the pion to be at rest {q = 0) in the medium, when it simplifies to 

d^p e{p,- - m 



n(")(,o,^-'=0) = -85i^go^|^ 



where n is the nucleon number density in symmetric nuclear matter, 



f2Ep -AEl + ql 
^^"^ ^o^ (3.13) 



Pf being the Fermi momentum, pp = \J 1^^ — fnif- Assuming that the vacuum part has already been taken care of to 

define the physical parameters, it is 11^"' (50) which is relevant in Eq. (3.7). 

It is now simple to calculate the remaining self-energy and vertex diagrams. A special feature is presented, however, 
by the constant vertex diagrams (f, g) arising from the vertex proportional to C2: While all other contributions are 
proportional to qxqa, this one is not, being given by 

(5"' iA{q){qxqaSaQ + g^go^Ao) 2c2n . (3.15) 

It reflects the fact that our treatment breaks Lorentz invariance to 0(3). Thus while the matrix element 

{0\A;\7T\q)) = l6'^"f,{q) , (3.16) 

is defined in vacuum as = (/^i^Tr, it must be expressed in medium as [23,9] 

U = <5^o qoFl + 6^, q^F^, z = 1, 2, 3 . (3.17) 

The results of calculating all the diagrams of Fig. 1 with vertices given by (2.17) may now be expressed in terms of 
the effective parameters, 

™<"l = m,{l+(2«-c.-c, + ^)^} , (3.18) 
i': = n{l+(« + c,-^)^}, (3.19) 

ir,. = n{l+(-c.+c,-^)^}. (3.20) 

These results were obtained earlier in this form in Ref. [10] by integrating out the nucleon field in the generating 
functional in presence of the external field a^. We postpone discussing the validity of these results until Sec. V. 
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IV. VIRIAL EXPANSION 



We next turn to a different approach to the problem, namely the virial expansion for the self-energy of the particle 
in question [4,11,12]. The resulting (first order) formula is valid, if the mcdimn is sufficiently dilute. As we shall 
discuss below, its range of validity is, in general, different from that calculated above using xPT. 

Let us derive the formula for the case at hand. A simple derivation follows, if we recognize that the self-energy 
function is an S-matrix element [12]. Consider first the process in vacuum. Just as the amplitude T for the two body 
ttN scattering, 

7r(fc, i) + N{p, s) 7r(/c', i') + N{p', s'), 

is given by the 5-matrix element, 

i{2T,f5\p + k-p'- k')Ti,i , = {k', i';p', s'\S - l\k, i;p, s) 

= {0\a{k', i'W, s'){S - l)6t(p, s)a^{k, i)|0) , (4.1) 

we may express the self-energy 11^°^ of the pion (in vacuum) by the one-body matrix element, 

-i{2T:f5^{k - k')5^,^\l^''\k) = {k',i'\S-l\k,i) = {0\aik',t')iS - l)a^(fc,i)|0) , (4.2) 
where S is the familiar scattering matrix operator, 

Here and s(s') arc indices denoting the pion isospin and the nuclcon spin projection in the initial (final) state 
respectively. Following the usual practice, the amplitude T is regarded as a 2 x 2 matrix in the nucleon isospin 
space. The operators o(fc, i) and b{p, s) annihilate respectively a pion of momentum k and isospin i and a nucleon of 

momentum p and spin s. 

The corresponding self-energy n(A;) in nuclear medium is obtained simply by replacing the vacuum expectation 
value in Eq. (4.2) by the ensemble average defined in Eq. (3.1), 

-i(2nfS^(k-k')Si,in{k) = {a{k',i'){S -l)aHk,i)) (4.3) 

It is here that we make use of the virial expansion in powers of the distribution function. We expand the ensemble 
average of any operator O as, 

{O} = (0|O|0) + ^ I ^^-||^n+(p)(p, s\0\p, s) + ..-, (4.4) 

where the sum is over the nucleon spin and isospin states. Applying this expansion to the right hand side of Eq. (4.3), 
we get 

-i{27rr6\k-k')IL^-\k)=Y^ I -^-^-^n+{p){p,s\a{k',i){S-l)a\k,i)\p,s) (4.5) 

where n(")(A;) stands as before for the difference, n(A;) — n(°)(A;). Notice that there is no sum over the pion isospin 
index i. We now use Eq. (4.1) to express the self-energy in terms of the forward scattering amplitude Tf{p, k), 



taking the summation inside the integral, as the distribution function is the same for all the four nucleon states. 

Eq. (4.6) is the desired first order virial expansion formula for the pion self-energy. Although such formulae have 
been used to find the mass shifts in different cases, its application to pion in nuclear medium does not exist in the 
literature. What has been utilised earlier is an approximation to the above formula [14] 
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But if the amplitude is not constant, even approximately, within the interval of integration in Eq.(4.6), this formula 
cannot clearly be trusted. As we shall sec below, this is indeed the situation for the case at hand. 

To sum over the forward amplitudes, we begin with the kinematics for the process in general. The matrix structure 
of T in nucleon isospin space may be written as 

TiH ,s's= TXSi'i + T^T, ^ , Ti] , (4.8) 

where each of the T"^ has the invariant spin decomposition, 

= u{p\ s'){A^ + m^Mp, s) . (4.9) 

Taking the amplitude for any one of the charged states n^'^ for the pion, we can readily carry out the sum over the 
nucleon states, 

Tf{p, k) = 8(mjvA+ + p ■ kB+) , (4.10) 

JV 

in terms of the isospin even amplitude only. 

We now evaluate 11'^"^ using the experimental data for the ttTV scattering amplitude [22]. As given by Eq.(4.6), it is 
an integral over the 3- momentum of the nucleon in the pion rest frame (fc = 0), while the data is given as a function 

of the pion energy = \JW^-rn^ in the nucleon rest (lab) frame. The two variables are related by the equation, 
m^Ep = mNOJk- We thus have finally the complex pole position in the pion propagator as 

_ = + s;^ = _ J, (^)' ^,,fj[z-.D*M ,4.11) 

where Z)+(a;fc) = A+ is the isospin even forward ttA'' scattering amplitude. The upper limit Wk is determined 

by the nucleon number density 

The imaginary part of the pole position represents the damping rate of pionic excitations in the medium. 

The numerical evaluations are shown in Figs. 2 and 3, where the pion mass shift and its imaginary part are plotted 
as a function of the nucleon number density in units of the normal density no=(110 MeV)^. Our result for the pion 
mass shift in nuclear medium may be compared with that for the nucleon mass shift in pionic medium, calculated 
in [4]. It will be observed that while both the shifts are given essentially by integrals over the same -nN amplitude 
times the corresponding distribution functions, the curves bend in opposite directions. The reason is that as the pion 
energy increases, there is a change in sign in the real part of the amplitude, which is weighted differently by the 
distribution functions in the two cases. 
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FIG. 2. Shift in pion mass in nuclear medium. The solid curve results from the virial formula. The three straight lines follow 
from : the dotted, the dashed and the dash-dotted curves correspond to the central value, the upper and the lower limits 
of the constants in Eq.(3.18). 
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FIG. 3. Damping rate of pionic excitations in nuclear medium. 
V. DISCUSSION 

Having derived the effective parameters of the pion in nuclear medium by two different methods, we now discuss the 
vahdity of the results. Considering xPT, the region in which the leading correction term may represent a meaningful 
approximation to a quantity depends on the proximity of the resonances in the relevant channel. In the present case, a 
number of resonances, particularly the A(1232), lie close to the threshold, making the values of the coupling constants 
ci , C2 and cs in the effective Lagrangian rather large. Then the region of nuclear density in which the results may be 
valid, is expected to shrink considerably. 

The calculated results follow this expectation. In the expression (3.18) for the effective mass, there is, however, 
a large cancellation among the contributions of the different vertices, making the mass shift to be only a few MeV 
at n = no, the normal nuclear density. But in the expression (3.20) for the contributions of the vertices add up, 
making the "correction" at this density overwhelm the unit term. Clearly, the first order results for the pion 

traversing nuclear matter at normal density are unacceptable. 

On the other hand, in the virial expansion formula we may avoid any inaccuracy in calculating the scattering 
amplitude by taking it from experiment. However, the nuclear medium must be dilute enough for an expansion of 
the self-energy function in powers of nuclear density to be valid. At normal nuclear density no = (110 MeV)^, the 
mean distance between the nucleons is about 2 fm. Thus our first order virial formula, where the pion propagation is 
perturbed by a single interaction with one of the nucleons in the medium, should be a reasonable approximation up 
to about this density. We point out that we use the exact formula, rather than the approximate version used so long 
by different authors in this context. Unfortunately, the virial expansion does not say anything about the residues, 
and F^ . At normal nuclear density, its prediction for the mass-shift and the decay width of the pion is that both are 
negligible, being only a few MeV. 

We conclude with a comment on determining the mass shift in the more interesting case of the nucleon in nuclear 
medium. Here the simple chiral Lagrangian for the NN system [25] is not expected to apply, as it does not take 
properly into account the presence of bound or virtual two nucleon states close to the NN threshold. In fact, it 
produces an absurdly large value for the nucleon mass shift [26]. But the virial expansion formula for the nucleon 
self-energy should apply here, at least in the neighbourhood of the nuclear saturation density [27] . 
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